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Abstract 

In this work we study certain invariant measures that can be associated to the time averaged observa- 
tion of a broad class of dissipative semigroups via the notion of a generalized Banach limit. Consider an 
arbitrary complete separable metric space X which is acted on by any continuous semigroup {S{t)}t>Q. 
Suppose that {S{t)}t>o possesses a global attractor A. We show that, for any generalized Banach limit 
LIM and any probability distribution of initial conditions tno, that there exists an invariant probability 

T— >CX3 

measure m, whose support is contained in A, such that 



(p{x)dm{x) = LIM^ f I (p{S{t)x)dmo{x)dt, 

T->oo7 Jo 



for all observables Lp living in a suitable function space of continuous mappings on X. 

This work is based on the framework of [FMRTOl]; it generalizes and simplifies the proofs of more 
recent works [Wan09],[LRRll]. In particular our results rely on the novel use of a general but elementary 
topological observation, valid in any metric space, which concerns the growth of continuous functions in 
the neighborhood of compact sets. In the case when {S{t)}t>o does not possess a compact absorbing set, 
this lemma allows us to sidestep the use of weak compactness arguments which require the imposition of 
cumbersome weak continuity conditions and thus restricts the phase space X to the case of a reflexive 
Banach space. 

Two examples of concrete dynamical systems where the semigroup is known to be non-compact 
are examined in detail. We first consider the Navier-Stokes equations with memory in the diffusion 
terms. This is the so called Jeffery's model which describes certain classes of viscoelastic fiuids. We 
then consider a family of neutral delay differential equations, that is equations with delays in the time 
derivative terms. These systems may arise in the study of wave propagation problems coming from 
certain first order hyperbolic partial differential equations; for example for the study of line transmission 
problems. For the second example the phase space is X = C([--r, 0], R"), for some delay r > 0, so that 
X is not reflexive in this case. 

Keywords: Dissipative dynamical systems, invariant measures, ergodicity, generalized Banach limits, global attractors, non- 
compact semigroups, systems with memory, Navier-Stokes, neutral delay differential equations. 
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1 Introduction 



Essential characteristics of dynamical systems are described by invariant probability measures. These mea- 
sures identify statistical equilibria and can provide important information about the long time behavior of 
the dynamics. It is therefore of paramount interest both in theory and applications to better understand 
this class of probability measures. 

Frequently however, statistically steady states represented by invariant measures are difficult to deter- 
mine. In practice physicists and engineers compute approximations of invariant probability measures by 
taking averages of time series (observables) associated to the system and invoking an 'ergodicity assump- 
tion'; they posit an equivalence between these temporal averages and averages against the unknown invariant 
measure which they are trying to identify. Mathematically speaking, the complete and rigorous justifica- 
tion of such an 'ergodic hypothesis' seems unreachable for non-specific classes of dynamical systems and 
remains a challenging problem even for specific examples. With this backdrop in mind, our aim in this and 
future works ([CGHb, CGHa]) is to further the development of a particular mathematical framework coming 
from [FMRTOl], which establishes a weak link between ensemble averages and temporal averages based on 
generalized limits but which we show is in fact applicable to a wide class of dissipative dynamical systems.^ 

For the study of infinite dimensional evolution equations a number of different mathematical approaches 
to ergodicity have been developed, each relevant to different situations. One approach has been to focus 
on some classes of simple linear or semi-linear first order partial differential equations. See, for example, 
[Las79, Daw83, BK84, Rud88, Rud04] and references therein. While these works are significant for the fact 
that they provide examples of linear systems where solutions exhibit chaotic or 'turbulent' behavior, the 
methods developed are limited to a very specific class of equations. Another, different, approach involves the 
study of stochastic partial differential equations (and other related infinite dimensional stochastic systems) 
where ergodicity is defined in terms of the Markov semigroup generated by the stochastic semifiow. See e.g. 
[DPZ96, FM95, KSOl, DPD03, HM06, HMll, Debll]. The methods developed in these works apply to a wide 
variety of stochastically perturbed dissipative linear and nonlinear systems, for instance the Navier-Stokes 
equations, Reaction-Diffusion equations. Delay Equations, weakly damped nonlinear Schrodinger equations, 
complex Ginzburg-Landau equations. On the other hand, this approach relics in an essential way on an 
underlying mechanism of stochasticity, one in which the noise injected into the system has to take a very 
specific form. Furthermore, this approach identifies invariant measures of the Markov semigroup which are 
deterministic probability measures and thus are not carried by the global random attractor.^ 

In this work we follow a different approach, pioneered in [FT75] and developed in [BCFM95; FMRTOl]. 
These works link ensemble and temporal averages via the notion of the so called 'generalized Banach limit', a 
linear functional acting over the space of bounded continuous functions which in particular associates those 
elements converging at -|-oo with their classical limit; cf. Definition 2.1 below. While the ideas in these works 
were developed in the specific setting of the 2D Navier-Stokes equations on a bounded domain some essential 
aspects of the framework generalize in a straightforward way to any compact semigroup {S{t)}t>o acting 
on a complete, separable metric space (X^d). On the other hand, in the important case of non-compact 
semigroups (see below for an extended discussion of examples of such systems) certain difficulties appear in 
generalizing this approach since the construction, as described in [FT75, BCFM95, FMRTOl], seems to rely 
in an essential way on the existence of a compact absorbing set. 

More recent works, [Wan09, LRRll], have been able to remove this compactness assumption by restricting 
X to be a reflexive Banach space and by requiring an additional 'weak-to-weak' continuity assumption on 
{5(t)}(>o. In essence, these works make use of the fact that bounded sets are weakly compact so that the 
arguments employed rely on the structure of reflexive Banach spaces and their associated weak topologies. 
While these requirements may be met in certain examples, the theory necessitates the verification of this 
additional weak-to-weak assumption which can be involved in practice (see e.g. [Ros98]). Moreover, there 

^ In this article, dissipative dynamical systems are defined as those which possess a bounded absorbing set, cf. [Tem97]. 
Note that other authors, e.g. [Hal88], have referred to dynamical systems with this property as having bounded dissipation. 

^On other hand these invariant measures are linked to the global random attractor; see [CSGll] for a discussion of such 
relationships in the context of stochastic differential equations (in finite dimensions). 
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are other interesting classes of dissipative dynamical systems which do not faU into this category. We develop 
some examples in detail below in this connection. 

We show here that these cumbersome assumptions, imposed in [Wan09, LRRll], are in fact totally un- 
necessary. We demonstrate that the corresponding results in [FMRTOl] extend to any continuous semigroup 
{S'(t)}t>o, evolving on any complete, separable metric space, which possesses a global attractor. Moreover, 
our method of proof is more elementary in character; we extend the methods of [FMRTOl] via a simple 
topological observation which limits the growth of continuous functions in a neighborhood of any compact 
subset of a metric space. See Lemma 3.1 below. 

With this background in mind we may describe the main abstract results in this work as follows: Suppose 
that {X,d) is any complete, separable metric space which is acted on by a continuous semigroup {S{t)}t>o 
possessing a global attractor A. We show that for any probability distribution mo of initial conditions over 
X, that there is an invariant probability measure m for {S(t)}t>o, which has its support contained in the 
global attractor A, such that 



for all real valued continuous bounded observables (p. Note that, as in [LRRll], we show here that (1.1) 
holds for an arbitrary probability distribution of initial conditions and not just for individual trajectories 
emanating from some xo & X; but again we do so without any restrictions on X beyond that it be a complete 
and separable metric space. On the other hand, in the case when mo is a Dirac measure supported on some 
point xo & X (i.e. when mo = S^o) wc show that (1.1) holds for any continuous real valued observable 
if. In any case, as with the previous results in this direction, (1.1) may be seen as establishing a kind of 
weak notion of ergodicity; by replacing the operation limr-j-oo with the operation LIM wc show that the 



equality of time averages and ensemble averages can be obtained. Note finally that our results generalize 
naturally to a non-autonomous or a stochastic setting but this comes at the cost of significant additional 
technicalities. We refer the interested reader to [CGHb, CGHa] where these generalizations will be carried 
out in the non-autonomous and stochastic cases respectively. 

The initial motivations that led us to discover the results appearing herein arose in ongoing work on 
nonlinear partial differential equations of parabolic type with memory effects added in the diffusion terms 
(see [CDPGHPll, CGHb, CGHa]). Such systems have a 'hyperbolic character' in comparison to their more 
classical cousins; they provide an interesting class of examples which do not possess in general a compact 
absorbing set for the associated semigroup. Of course there are many other important evolution equations 
which exhibit non-compact semigroups; for example, the Navier-Stokes equations on unbounded (i.e. domains 
where the Poincare inequality holds) [Ros98] , retarded equations with infinite delay [HL93] , neutral functional 
differential equations [HL93], certain partial functional differential equations [Wu96], the linearly damped 
nonlinear wave equation as well as other equations of (partially) hyperbolic type [Lad91, Tcm97]. Note that 
restricting consideration to semigroups evolving on reflexive Banach spaces is too stringent for many of the 
above cited equations. For instance in the case of neutral functional differential equations, cf. [HL93], the 
phase-space is typically C([— r, 0],M") for some delay time r > 0. 

Here, in order to illustrate the flexibility of our main abstract results, we will study two dynamical systems 
in detail. In each case we consider systems with memory which generate non-compact semigroups. We first 
consider a model for a viscoelastic fluids similar to the 2D Navier-Stokes equations but incorporating non- 
local, integro-differential diffusive terms that depend on the past history of the flow. See e.g. [Jos90, AS98, 
Orl99, GGP05] and below for further mathematical and physical background. We then turn to a class of 
neutral delay differential equations (NDDEs) . These are functional differential equations with dependence on 
the past values of the solution and its time derivatives. Such systems arose in relation to the study of certain 
line transmission problems modeled by the telegrapher's equation [BM64], and have been since encountered 
in various engineering and physical applications involving other hyperbolic PDEs; see e.g. [KHIO] and 
Remark 4.5 below for connections to systems in geophysical fluid dynamics. For further background on the 
mathematical theory of NDDEs see e.g. [HL93]. 




(1.1) 



T-^oo 
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The manuscript is organized as follows: We begin by briefly reviewing some mathematical generalities 
and setting notations that will be employed throughout the rest of the work. The main results are then given 
in precise terms in Theorems 2.1 and 2.2 below. Complete, self-contained, proofs of both results are next 
given in Section 3. Finally in Section 4, we turn to the study of the two concrete examples, the Navier-Stokes 
equations with memory and certain classes of NDDEs, establishing novel results linking invariant measures 
to temporal averages for these systems in Theorem 4.1 and Theorems 4.2, 4.3 respectively. 

2 Notations, mathematical preliminaries and statement of the 
main results 

Before stating the main abstract results of the work in precise terms (Theorems 2.1, 2.2) we first review some 
basic definitions and other essential mathematical preliminaries setting notations that will be used below. 
Throughout the rest of the article we will always take {X, d) to be an arbitrary complete, separable metric 
space and consider a continuous semigroup {S(t)}t>Q on X] more precisely we assume that 5(0) = Idx, 
S{t + s) = S{t)S{s) for all t,s £ M+, and that S : x X ^ X is separately continuous. 

Recall that a global attractor ^ is a compact subset of X that is invariant under S{t), i.e. such that 
S{t)A = A for alH > and which attracts all bounded subsets of X viz. 

lim dHiSit)B,A) = 0, for all B ClX,B bounded, 

where dn is the HausdorfF semi-distance 

dniE, F) := sup inf d{x, y) for any E,F C X. 

The study of attractors is an extensive and well-developed subject see e.g. [Tem97, Hal88, Lad91, Chu99, 
RobOl, CV02, MZ08]. For the abstract results. Theorems 2.1, 2.2, we will assume that {S'(i)}t>o possesses 
a global attractor A. On the other hand, for the concrete examples considered in Section 4, we employ the 
following useful and rather general sufficient condition for the existence of a global attractor. Note also that 
this result follows immediately as a special case of the characterizations of semigroups possessing a global 
attractor appearing in e.g. [MWZ02, Theorem 3.8] or [CCP, Theorem 11]. 

Proposition 2.1. Let H he a Banach Space (with an associated norm \\ ■ \\) and consider {S'(t)}t>o a 
continuous semigroup acting on H . Suppose that 

(i) there exists a hounded set *B C -ff such that for every B d H , B bounded, there exists t* ~ t*[B) > 
such that S{t)B C 05 for every t > t* . 

(ii) For t > we may split S{t) as S{t) ~ Si{t) + S2{t) such that, for every K > 0, 

sup ||5i(i)a;|| 

xeH:\\x\\<K 

and for every bounded set B and every t > 0, S2{t)B is a precompact subset of H . 
Then {S{t)}t>o has a (connected) global attractor A which is the omega limit set of % i.e. 



A = u:m ■.= {^[jS{s)^. 

t>0 s>t 

We next introduce some notations and recall some basic properties associated with probability measures 
defined on the general metric space X. Let Pr{X) be the set of all Borel probability measures on X with 
B{X) the associated collection of Borel measurable sets. For m G Pr{X), we define supp(m) to be the 
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smallest closed set E such that m{E) = 1. See e.g. [Rud87].'^ Take C{X) (resp. Cb{X)) to be the collection 
of real-valued contmuous (resp. bounded continuous) functions defined on X . A measure m G Pr{X) is said 
to be invariant (relative to {S'(i)}t>o) if 

m{E) = m{S{t)~^E), for alH > and every E e B{X), (2.1) 

or equivalcntly if 

/ ip{x)dm{x) = / v5(S'(<)x)dm(a;), for alH > and every e C(,(X). (2.2) 
Jx Jx 

Note that, for every invariant measure m, supp{m) is contained in the global attractor A. For completeness 
we recall the proof of this fact, which is elementary, in an Appendix. See Lemma 4.7 below. 

Recall that a sequence {m„}„>o C Pr{X) is said to converge weakly to a measure m ifFlim„_j.oo ^dmn = 
Jj^ ipdm, for every (p G Cb{X). As such {Tn„}„>o may be said to be weakly compact if we can extract from 
{iTi,i}„>o a weakly convergent subsequence. On the other hand we say such a collection {m„}„>o is tight 
if, for every e > there is a corresponding compact set C X so that mn{K^) > 1 — e, for every n. 
Classically these two notions, tightness and weak compactness, are equivalent, a result usually referred to as 
Prokhorov's theorem. See e.g. [Bil99]. 

Let us also recall a special case of the classical Kakutani-Riesz Representation theorem, as suits for our 
purposes below. See e.g. [Rud87] for further details. 

Lemma 2.1. Let K he a compact Hausdorff space and suppose that £, is a positive linear functional on C{K) 
(the continuous real valued functions on K with the usual sup norm). Then there exists a unique positive 
Borel measure m on K such that, for every (p G C{K), 2,{ip) = ip{x)dm{x). 

Finally we turn to the notion of a generalized Banach limit which is defined as follows: 

Definition 2.1. Consider the collection B+ of all hounded real-valued functions on [0,oo) endowed with sup 
norm. A generalized Banach limit, which we denote hy LIM, is any linear functional on S+ such that 

(a) hlMgU) > for all g e B+ with g{s) > 0, for all s > 0. 

T— >oo 

(b) lAMg{t) ~ lim g(t) for all g G for which the usual limit exists. 

T^oo t^oc 

It is not hard to establish the existence of such a positive linear functional as a consequence of the Hahn- 
Banach theorem. Note also that, for any such LIM, it may be shown that 

T-s-oo 

LIMg(r) <limsup|5(T)|, (2.3) 



T->oo 



for any g G B+. Wc will use this observation frequently below. Sec e.g. [Lax02] for further background and 
properties. 

With these preliminaries in hand wc now state the first main result which shows that (1.1) holds for any 
continuous observable in the case when mo is a Dirac measure. 

Theorem 2.1. Suppose that {X,d) is a complete, separable metric space and {>S'(t)}f>o is a continuous 
semigroup on X that possesses a global attractor A. Fix a generalized Banach limit LIM . Then, for any 

T~yoo 

xo G X , there exists a unique invariant measure m G Pr{X) for {S{t)}t>o such that 

p{x)dm{x) = UM^ [ p{S{t)xo)dt, for any p € C{X), (2.4) 
and such that supp{m) C A. 



X T^ooT jQ 



^ Note that a related, but purely measure theoretic notion of a 'carrier' is also sometimes used in this connection. However, 
both as it concerns us here and in its usage in previous related works, cf. [FMRTOl, BCFM95], the two notions can be 
seen to be applied in a completely equivalent fashion. We will make this equivalence precise after the statement of the main 
Theorems 2.1,2.2 in Remark 2.1, (iii) below. 



5 



The second result establishes (1.1) for any initial probability measure mo but requires further in this generality 
that the observable be both continuous and bounded. 

Theorem 2.2. Suppose that {X^d) is a complete, separable metric space and {S{ty\t>o is o. continuous 
semigroup on X that possesses a global attractor A. Fix a generalized Banach Limit LIM . Then, for any 

T— foo 

TTio € Pr(X), there exists a unique invariant measure m G Pr{X) for {S{t)}t>o such that 

[ ip{x)dm{x) = LIM- [ I Lp{S{t)x)dmo{x)dt, for any Lp e Cb{X), (2.5) 

Jx T^ooT Jq Jx 

and such that supp{m) C A. 

We conclude this section with some remarks concerning Theorems 2.1, 2.2. 
Remark 2.1. 

(i) The uniqueness of the invariant measures m satisfying (2.4) or (2.5) follows as a direct consequence of 
the regularity of borel probability measures on metric spaces; see e.g. [Bil99]. We therefore need only 
to establish the existence of such measures in the proofs below. 

(a) By considering, for xq € X , the Dirac measure xviq = S^o we partially recover Proposition 2.1 from 
Proposition 2.2. Notwithstanding, we separate the two results since, in the former case of Proposi- 
tion 2.1, we are able to establish (2.4) in the larger class of test functions C{X). Note that the use 
of the smaller collection Cb{X) for the space of test functions in Proposition 2.2 allows us in partic- 
ular to work in the topology of weak convergence of measures. This is needed to be able to pass to a 
limit in a sequence of approximating measures via the Prokhorov Theorem. See Step 2 of the proof of 
Proposition 2.2 below. 

(Hi) Theorem 2.1, in the given generality of complete separable metric spaces acted on by arbitrary semi- 
groups possessing a global attractor, also appears in [LRRll]. We give a different proof of Theorem 2.1 
below, based on Lemma 3.1. Lemma 3.1 is also used in an essential way in the proof of Theorem 2.2. 
In contrast to Theorem 2.1, Theorem 2.2 establishes results analogous to those appearing in previous 
works in a much greater generality and is thus new. 

(iv) In other related works, e.g. [FMRTOl, BCFM95], the notion of a carrier is sometimes used as an 
alternative to the support of a measure. Recall that an element m € Pr{X) is said to be carried by 
a set E G B{X) (or that E is a carrier for m.) if m{E) ~ I. Of course, m is carried by supp(m) (the 
smallest closed set of full measure ), but this particular carrier is not unique in general ( and rarely 
would be in practice). In any case, since the global attractor A is closed, when say in Theorems 2.1, 
2.2 that 'supp{m) ^ A' we may equivalently state that 'm is carried by A'. 

(v) The so called Krylov-Bogoliubov procedure provides another means of associating invariant measures 
with time averages starting from a given fixed initial measure fiQ. In [LRRll] it was shown that this 
procedure could be used to establish the existence of invariant measures for the class of dynamical 
systems considered in, for example. Theorems 2.1, 2.2. 

3 Proof of abstract results 

We turn in this section to the proof of Theorem 2.1, 2.2. Both proofs lean heavily on the following general 
topological lemma which is established along elementary lines. 

Lemma 3.1. Let {X,d) be a metric space and consider any compact subset K C X. For e > Q we let, 
Ke :~ {x ^ X : inf j,gA' d{x, y) < e}. Then the following properties hold: 
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(a) For every (p G C{X), there exists e > so that sup < oo. 

(h) Suppose ifjip € C{X) are such that (p{x) = ^pi^) for every x G K, then for every e > there exists 
i5 > such that, sup \f{x) — ip{x)\ < e. 

x^Ks 

Proof. For the first item, (a), fix (^9 G C{X). For every x £ K wc may choose S — 6x > so that, for 
every y G B{x,Sx) := {y : d{x,y) < Sx}, I'fiix) — ^{y)\ < 1- Choosing numbers Sx > in this manner 
we may form the open cover C = {i?(.T, (5j;/3) : x G K\ for K. Since K is compact, we can extract 

from this cover a finite sub-cover C = {i3(a;i, (5^:^/3), i?(a;„, (5^:^/3)}. Take e = """^'^^^g'"''^''"'^ and let 
i\/ = f + maxj=i^...^„ \ip{x.j)\. 

Given any x G we may choose y & K such that d{x^ y) < 2e. Since C covers K we may pick Xj such that 
d{y,Xj) < 5xji- Combining these two observations we conclude perforce that d{x,Xj) < 2e + 6x^/5 < Sxj 
so that |¥'(a::)| < M. Since x £ was arbitrary to begin with this gives (a). 

We turn to the second item (b). Fix e > 0. For every x G -K' we choose > so that \<f{x) — f{y) \ + 
\^l^{x) — V'(y)| < e whenever y G B{x,jx)- Again, due to the compactness of K, we may cover K with a 
finite collection {B{xi,jxi /3),...,B(x„,7^„/3)}. Take S = '""'^^"^g"-^""^ Similarly to the previous case 
we observe that Kg C ^k=i-^i-^k,7k)- Fix arbitrary y G Kg and choose k so that y G B{xk,^k)- Noting that 
(p{xk) = ij{xk): this implies \(p{y)-il;{y)\ = \ip{y)-(p{xi) + ip{x^)~tp{y)\ < \Lp{y)-ip{xi)\ + \tl;{xi)~tp(y)\ < e, 
as needed for (b). The proof is therefore complete. □ 



Proof of Theorem 2.1: 

Fix any xq € X . We proceed in steps. First we show that the generalized Banach limit on the right hand 
side of (2.4) exists. We then show that the resulting functional is uniquely defined by the restriction of ip to 
the global attractor A. Since A is compact, this is sufficient to infer the needed m via an application of the 
Kakutani-Riesz theorem. In the final step we show that the m that we have found is indeed invariant. 
Step 1: Existence of the positive linear functional. To carry out this first step, wc show that the map 
defined by T H> y ip{S{t)xQ)dt is bounded over [0, 00) so that 



£,„(<^) LImI / ip{S{t)xQ)dt, for(^GC(X), (3.1) 



is well defined as a positive linear functional on C{X) (cf. Definition 2.1). Since, by assumption {S'(t)}t>o 
possesses a global attractor A, we have that, 

for every e > 0, there exists a time > such that S{t)xo G A^ for every t > T^, (3.2) 

where here and below, Ae = {x € X : inf d(x, y) < e}. Applying Lemma 3.1, (a) to A and (p we 
may choose e > so that Ka sup2.g_4^ l¥'(a;)| < 00. Taking as required in (3.2) for this value of e 
we find also that Kj := supjgjoT,] l'P(*5'(i)^o)| < 00 owing to the fact that the interval [0,rf] is compact 
and that t i~>- S{t) is continuous. By now taking K = maxjif/, Xa}, we infer, for every T > 0, that 
T lo ^{S{t)xo)dt < K < 00 as needed to finally justify the definition of 2,xa given in (3.1). 

Step 2: Restriction to A and the application of the Kakutani-Riesz theorem. The next step will be 
to show that (ip) depends only on the values of (p on A. More precisely, we establish that if (p{x) = (p{x) 
for every x £ A then £2:^ ((p) = &xo (<^). To this end we fix any e > and according to Lemma 3.1, (b) choose 

a corresponding 5 > such that sup Ifix) — ip{x)\ < e. We now take Ts so that S{t)xo G As for each t > Ts 

xeAs 

and let Kg '■= suptgjg rp^^{\ip{S(t)xo)\ + \'p{S{t)xo)\). For the reasons noted for Ki in the previous step above. 
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Ks is finite. Using the basic properties of LIM, (2.3), we estimate 



LIM- ^{S{t)xo) ~ ^{S{t)xQ)dt 







< linisup 







ip{S{t)xo) - ip{S{t)xo)dt 



< lim sup — 

T->-oo T 





< lim sup — — — 

T-S.OO T 



\(p{S{t)xo) ~ 'f{S{t)xn)\dt + lim sup - 

T^oo J 

T 



\ip{S{t)xo)-ipiSit)xo)\dt 



Ts 



1 



limsup- / \ipiSit)xo)-^{S{t)xa)\dt 



< lim sup h iini sup < e. 

T-i-oo 2^ T^cxj r 

Since, to begin with, the choice of e > was arbitrary wc conclude perforce that £>xo{v — = as desired. 
With this in hand we may now unambiguously define 

©(V)=£xoWV')), iori,eC{A), (3.3) 

where we take £ : C{A) — >■ C{X) to be an extension operator, £{il}){x) — ip{x) for x & A, such that 

inf £{ip){x) = inf (p{x), sup £{(p){x) = swp ip{x). (3-4) 

The existence of such an extension operator is guaranteed by the Dugundji extension theorem, [DugSl, 
Theorem 4.1]. One may readily verify that © is linear. With (3.4), (3.1) and Definition 2.1, (a) we see that 
<& is also positive. Thus, by the Kakutani-Riesz representation theorem, recalled above in Lemma 2.1, there 
exists a unique positive, finite, Borel measure m on ^ such that 



©(V) 



ip{x)dm{x), for ip G C{A). 



(3.5) 



Abusing notation slightly, we extend m to a Borel measure on all of X by taking m{E) := m(£' n A), 
E e B{X). Clearly, m(\4=) = 0, and so with (3.1), (3.3), (3.5) it follows that, for every tp e C{X), 
Zxoif) = /_4 'p{x)dm{x) = j-^ Lp{x)dm{x) , which is (2.4). To see that m lies in Pr{X) we simply insert 1,5 = 1 
into (2.4). 

Step 3. The invariance of m. For the final step we show that the m we have found in Step 2 is indeed 
invariant by establishing (2.2). To this end fix any t* > and any E Cb{X). We compute: 



(p{S{t*)x)dm{x) = LlM^ [ ip{S{t*)S{t)xo)dt ^ LIM^ [ ip{S{t + t*)x)dt 

X T~>ooT Jq T^ooT Jq 

T+t* 

(p{S{t)xo)dt 



1 



= LIM- 



1 



T+t* 



(^(x)dm(x) + LIM- 
•p{x)dm{x). 



^(5(i)xo)dt- limI 

T— >oo 1 



t* 



(3.6) 



ip{S(t)xQ)dt 



X 



Note that the first equality is justified since i/'(') = € Cb{X). For the final equality above we note 

that both LIMtJt C ip(S(t)xo)dt and LIM ^ /tT^* 'fi(S(t)xo)dt are zero due to the boundedness of ip and 
the basic properties of LIM , Definition 2.1, (b). We have thus completed the proof of Theorem 2.1. 

T— >oo 
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Remark 3.1. For the proof of Theorem 2.1, the fact that supp{xn) C A follows directly from the application 
of the Kakutani-Riesz theorem on the global attractor A. The situation is different in e.g. [FMRTOl] where 
Kakutani-Riesz theorem is applied on a compact absorbing set. Here one needs to rely on the general result, 
recalled in Lemma ^.7 and cf. [FMRTOl, Chapter 4, Theorem 4-1], that the support of an invariant measure 
is always contained in the global attractor A. On the other hand we use Lemma ^.7 in Step 2 of the proof of 
Theorem 2.2 for the general case of an mo G Pr(X) with an unbounded support. 



Proof of Theorem 2.2: 

Let m-o e Pr{X) be given. Since we are assuming here that ip is bounded, in contrast to the proof of 
Proposition 2.1 it is trivial to define: 

£n,oM = LIM;J; / ( p{S{t)x)dmQ{x)At, iovpeCb{X). (3.7) 

As in [LRRll], we now proceed to establish (2.5) in two steps. Initially we assume that mo has a bounded 
support so that this support is attracted to A. At the second step we drop this assumption and pass to 
the general case of any mo G Pr(X) with compactness arguments involving the Prokhorov theorem and a 
suitable sequence of probability measures mo approximating mo. 

Step 1: mo with bounded support. Assuming that supp(mo) is a bounded subset of X, we fix any 

(p,if> ^ Cb{X) with (f{x) ~ (fix) for all x G A. As in Step 2 of the proof of Proposition 2.1 we would 
like to show that Slmoiy^) = ^moi'P) so that we may once again apply the Kakutani-Riesz theorem. For 
this purpose fix e > 0. Invoking, as above. Lemma 3.1, (b) we choose a corresponding 5 > such that 
sup2,g_4^ \ip{x) — ip{x)\ < e. Since A is attracting and supp(mo) is assumed to be bounded, we may choose 
then choose a time t* > such that 5'(i)(supp(mo)) C As for all t > t* . Using basic properties of LIM, 

(2.3), we estimate 

|£m„(</'-<^)| <limsup;J; / / \ip{Sit)x)-ip{Sit)x)\dmoix)dt 

T->-oo J Jo Jx 
1 '•^ 



+ limsup— / / \ip{S {t)x) - ip{S{t)x)\ dmo{x)dt 

T-i-oo T J ft J 

supp(mo) 

rsup,g^(|y>(x)| + |^(x)|) , ^. {T^t*)e 
< hm sup — h um sup — = e. 

Since e was arbitrary to begin with we infer that ZmtX^) = '2mo(<^) for all <p, <^ e Ch{X) with Lp\j( = ip\j(. 

With this in hand we are able to define a functional © as in (3.3) and obtain an associated m with 
m(./l'^) = 0.''^ With this measure m, (2.5) now follows exactly as in Step 1 of the proof of Proposition 2.1. To 
show that m G Pr{X), we take (p = 1 m (2.5) and use that mo G Pr{X). 

Step 2: passage to the general case via the Prokhorov theorem. Now suppose that mo G Pr{X) 
is arbitrary. Since mo is a Borel probability measure it is tight] we may choose a sequence of compact sets 
Kn C X such that mo(i^n) > 1 — 1/n (see e.g. [Bil99]). Accordingly we now define a sequence of measures 
m[)' Pr{X) via 

mo{Ef^Kn) 



Clearly each mo has a bounded support. Invoking Step 1, we obtain, for each n, a corresponding invariant 
measure m" such that mfj and m" satisfy (2.5). Note also that, for each n supp(m") C A. Hence, since 
A is by definition a compact set, by the Prokhorov theorem (see e.g. [Bil99]), we infer that {m"}„>i must 



*Notc that, due to (3.4), the extension operator £ is chosen to map C{A) into Ch{X). 



9 



have a weakly convergent subsequence. Thinning this sequence of m"'s as necessary we find that, for each 



f 11 

(p{x)dm{x) = lini / ip{x)dm"{x) = hm — -LIM— , , 

T 



ip{S{t)x)dmQ{x)dt 



¥ ;;nirT (^I^?^ / / ^(5W^)dmo(x)di- LIM^ / / ^iSit)x)dmo{x 



)dt 



LIM- I / (f{S{t)x)dmo{x)dt, 
10 Jx 



which is (2.5). Note that we may justify the last equality above by observing that 



ip{S{t)x)dmo{x)dt 



Jx\K^ 



< sup \(p{x)\mo{X \ Kn) < 
xex 



We show that m G Pr{X) exactly as in the previous step by inserting (p = 1 into (2.5). The invariance of 
m (for the case of a mo with bounded support or otherwise) follows from a computation like (3.6). Having 
shown that m is an invariant probability measure it follows that supp(m) C A. Indeed this is a general 
property of invariant measures that we restate and prove for the sake of completeness in Lemma 4.7 below. 
The proof of Theorem 2.2 is now complete. 

4 Application to evolution equations with memory 

In this section we consider two concrete dynamical systems which are non-compact but which nevertheless 
possess a global attractor. The common theme between of these examples is the presence of memory terms 
in the governing equations; we suppose that the evolution of the state variables depends on both past and 
current states the of system. As noted above in the introduction, each of these examples are not covered 
under the previous results appearing in [Wan09], [LRRll]. In the first section we study a variation on the 
Navier-Stokes equations, often referred to as the Jeffery's Model, which has diffusive memory terms. Further 
on in Section 4.2, we consider a class of neutral differential equations, i.e. differential equations involving 
time derivatives of the state variable at lagged times. Here, in addition to the fact that the semigroup is 
non-compact, the underlying phase space X is a non-reflexive Banach space. 



4.1 A viscoelastic fluids model: The Navier-Stokes equations with memory 

For the first example we consider a variation on the Navier-Stokes equations, the so called Jeffery's models 
which incorporates the past history of the flow through additional diffusive, integro-differential 'memory' 
terms. See (4.1) - (4.4) below for the precise formulation of this model. Such equations arise in the study 
of certain viscoelastic fluids; for example to model dilute solutions of polymers or bubbly liquids. See 
[Jos90, GGP05]. 

Physically speaking, viscoelastic materials exhibit effects of both elasticity and viscosity. For such ma- 
terials, the stress is typically a functional of the past history of the strain, instead of being a function of 
the present strain value (elastic) or of the present value of the time derivative of the strain (viscous). When 
an integral term is used for the history dependence, the dynamical equations become partial integrodiffer- 
ential equations with a character somewhere between hyperbolic (elastic) and parabolic (viscous); see e.g. 
[FS86, RHN87]. Thus, it is not surprising that the Jeffery's model has interesting hyperbolic properties in 
comparison to the classical Navier-Stokes equations, as we will see below. 

Numerous models coming from mathematical physics or mathematical ecology incorporate such memory 
terms. For example one may add memory to the damped or strongly damped wave equation for the study of 
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wave propagation in materials with memory, or to nonlinear parabolic equations which result in "reaction- 
diffusion systems" that take non-Fickian diffusive effects into account. Such parabolic equations are used in 
the modeling of heat propagation in certain materials [CG67, GP68] or more recently have been advocated 
in the study of population dynamics in situations where the individuals spread according to environmental 
feedbacks, see e.g. [Pao97, CG06]. 

In this type of "reaction-diffusion system" with non-Fickian diffusive effects, the contribution of the 
medium to the change in mobile concentration is modeled by a linear density exchange process which may 
be viewed as a source term with respect to the memory-less (i.e. Fickian) reaction-diffusion system. This 
source term can be expressed as the convolution product of a time-dependent function k, called the memory 
kernel, and some spatial variation of the mobile concentration u. For instance this might lead to a diffusive 
memory term of the form J_ K,(t — s)Au{s)ds. In any case, main feature of this formulation is that the 
memory kernel depends only on the properties of the medium, and is therefore an intrinsic characteristic 
of the system. See [GMLB+08] for other types of convolution products between the time derivative of the 
mobile concentration and the memory kernel as arising in non-Fickian dispersion in porous media and other 
'real-world' situations. 

The mathematical study of hyperbolic and parabolic systems modified by the kind of diffusive memory 
terms discussed above is now extensive. See e.g. [Daf70, GP06, CDPGHPll, PZOl, GMPZ08, GP05, 
GGMP06, DPPZ08, DPP08, CPS06, CP06, CGG+06, GGP99, GP02, GP]. The Jeffery's model we will 
consider here was previously studied in [AS98, Orl99, GGP05]. Since we are interested in the dynamical 
properties of this system, we will treat (4.1) - (4.4) as an autonomous dynamical system within the extended 
phase space formalism developed in [Daf70] and used subsequently in many of the works just cited. While 
the existence of an attractor for (4.1) - (4.4) in this extended phase space setting was studied in [GGP05] 
we will establish such results here under a larger class of conditions on the memory decay kernel k than was 
possible in [GGP05] . We are able to accomplish this feat by making use of an appropriate energy functional 
(see Lemma 4.2 and also [CDPGHPll, GMPZ08]). Note that while this more general condition, (4.6), is 
much more challenging technically, it allows for the interesting case of a memory term k exhibiting a linear 
decay as we describe below in Remark 4.1. Having established that (4.1) - (4.4) possesses a global attractor 
within this framework the existence of invariant measures naturally follow from Theorem 2.1 and Theorem 
2.2. In any case, to the best of our knowledge, no one has previously considered the existence of invariant 
measures for the Jeffery's model, or for that matter any other parabolic or hyperbolic systems with such 
diffusive memory terms. 

Note that the time asymptotic dynamics of (4.1)-(4.2) can be quantified in the cases of a time dependent 
or a stochastic, white noise type forcing. This will be carried out elsewhere in [CGHb, CGHa]. Of course, 
we could formally consider a three dimensional version of (4.1)- (4.2). However, as with the classical three 
dimensional Navier-Stokes equations, the existence and uniqueness of strong solutions is unknown. The three 
dimensional system is thus far from the reach of the theory we have developed here. 

4.1.1 The Jeffery's model: governing equations and basic assumptions 

We introduce the model as follows. Fix a bounded open domain 2? C with smooth boundary dV. On D 
we consider the following integro-differential equation 

dtu + u-'Vu^iyAu+ f k(< - s)Au(s)ds + Vp + /, (4.1) 
V-M = 0. (4.2) 

Here u = {ui,U2), p represent the flow field and the pressure of a viscous incompressible fluid filling V. 
The function k determines the dependence on the past history of the flow through the integro-differential 
'memory' term K{t — s)Au{s)ds and distinguishes (4.1)-(4.2) from the classical Navier-Stokes equations. 
As such, K is assumed to be positive and decreasing. Further decay conditions for k are be imposed below in 
(4.5)-(4.7). We will assume throughout what follows that the external body forcing / is time independent 
and that / G L'^(V). 
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Of course (4.1)-(4.2) is supplemented with initial and boundary conditions. We impose the no-slip 
(Dirichlet) boundary condition 

u\av = 0. (4.3) 

Observe that, in view of the history term n{t — s)Au{s)ds, the divergence-free vector field u must be 
known for alH < in order to make sense of (4.1). Accordingly, the boundary- value problem is supplemented 
with the 'initial condition', or more precisely the initial past history: 

uit) = uo(i), t < 0. (4.4) 

Of course for each t < 0, uo{t) respects the conditions (4.2), (4.3). 

Concerning the memory terms in (4.1), we assume that k may be written in the form: 

ix{(j) da := — I ii{a)da, (4-5) 

for some nonnegative, nonincreasing function ^ G L^(M+), so that n is nonincreasing and nonnegative. We 
will sometimes refer to this /z as the memory kernel as well, when no confusion is possible with k. Note that 
k(0) = kq is the 'total mass of ^' i.e. Jp°° ^(s) ds = kq and that k' = — /i for almost every t > 0. Throughout 
the following, we impose the decay condition on fi 

fj.{s + a) < Ke^^^nis), a.e. s,cr>0, (4.6) 

for any fixed K > 1, S > desired. This is equivalent to requiring that 

k{s) < Ptiis), a.e. s > 0, (4.7) 

where we may take [3 ~ As a byproduct of (4.7), k £ _L^(M+). Using additionally that —k' = wc have 

k(s) < Koe""/'^ = Koe^^'^^/^"^ for every s > 0. 

Remark 4.1. The condition (4.6) is more general in comparison to many previous works including [GGP05] 
where it was required that K = I. The 'physical' significance of (4.6) is that it allows the treatment of a 
memory kernel /i that has 'flat zones', i.e. (4.6) allows for the consideration of a k with an linear decay 
profile. Indeed, take 

, N J Mo for s < t*, 
' [0 fors > t*, 

where t* , no > are fixed constants. Then, according to (4.5), kq = ^ot* and 



K[S 



{l-s/t*)Ko fors<t* 
for s> t* 



Observe that (4.6) is satisfied by any d > by taking K = e^^ > 1. On the other hand (4.6) is clearly 
violated for K ~ 1, regardless of the choice of S > 0. 

4.1.2 The extended phase space and its functional setting 

In order to treat (4.1) - (4.4) as an autonomous dynamical system we follow [DafZO] (and sec above for 
extensive further references) and introduce an additional 'memory' variable rj, which is defined according to 



r]\s) := [ u{t - a)dcr = / u{a)da, for t > 0, s > 0. (4.8) 

Jo Jt-s 
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Arguing formally, we obtain from (4.1) - (4.4) the following coupled system of equations 

dtu + u-'^u = i^Au- fi{s)AT]*{s)ds + \/p + f (4.9a) 
Jo 

dtv*{s) = -dsv\s)+u{t), (4.9b) 

V • u = V • ry = 0, (4.9c) 

u{0) = uo{0), 7fis) = / ^0(^)01^, (4.9d) 



u\av = 0, -nidv = (4.9e) 

Wc recall this derivation in detail [CGHb]. 

In order to place (4.9a)-(4.9e) in a rigorous functional framework we next recall some classical spaces in 
the mathematical theory of the Navier-Stokes equations. See e.g. [TemOl] for further background. Further 
spaces needed for the memory variable 77 will be recalled further on below. 

Take U -.^ {(t) e {C^{V)f : V • = 0} and define H := dL^^v)^^ = {u e L^{Vf : V • u = 0, u • n = 0}. 
Here n is the outer pointing normal to dD. On H we take the inner product {u, v) := J-pU ■ vdD and 
associated norm \u\ := \J {u, u). The Leray-Hopf projector, Ph, is defined as the orthogonal projection of 
L'^iV)^ onto H. At the next order let V -.^ cIh^(v)U = {u G HI{'DY : V • u = 0}. On this set we use 
the inner product {{u,v)) := Vu • S/vdV and norm ||m|| {{u, u)). Note that due to the Dirichlet 

boundary condition the Poincarc inequality \u\ < Aj~"'^||u|j holds for all u € V. Here, the constant Ai is the 
first eigenvalue of the stokes operator A defined in the next paragraph. This justifies taking || • || as a norm 
for V. We take V to be the dual of V, relative to H with the pairing notated by (•, •). 

We next define the Stokes operator A which is understood as a bounded linear map from V to V' via: 

{Au, v) = {{u, v)) u,veV. (4.10) 

A can be extended to an unbounded operator from H to H according to Au = —PhAu with the domain 
D{A) = H'^{pY n V . By applying the theory of symmetric, compact, operators for A~^ , one can establish 
the existence of an orthonormal basis {ek\k>i for H of eigenfunctions of A. Here the associated eigenvalues 

{Afe}fe>i form an unbounded, increasing sequence viz. < Ai < A2 < . . . < A„ < A„+i < We shall also 

make use of the fractional powers of A. For u € H, wc denote Uk = (w, e^). Given m > take D{A™'/^) = 
{uGH:J2k ^ri"fcP < 00} and define A^/^w = J^k >^T^'^Ukek, u e DiA""). We equip D{A^''^) with the 

norm |w|,„ := = {Y.k K\uk?f'^ ■ Note that V D{A^'^) and that for u e D{A^'^), \\u\\ = \u\i. 

The nonlinear portion and part of the pressure gradient in (4.9) is captured in the bilinear form: 

B{u, v) := Ph(u ■V)v= ^ PniujO-jv) u e V, v e D{A). (4.11) 



For notational convenience we will sometimes write B{u) :— B{u,u). Note that B is also well defined as a 

f^{u ■ V)v ■ wdV = Ej,fc=i Iv '' 



continuous map from V ^ V to V according to {B{u, v), w) :— J-p{u ■ V)u • wdD — f.^^ UjdjVkWkdT). 



It is easy to show that 

(B(m, w),u) = 0, forallu,wel/ 

Classically, with Holder's inequality and Sobolev embeddings we have the estimate 

\{B{u,v),w)\ < C|u|i/2||y||i/2||^|||y^|i/2||^^,||i/2^ fQ^. u,v,we V. (4.12) 

On the other hand, slightly different estimates along the same lines lead to, 

|(B(w,u),ii;)| < C|ii|^/2||w||^/2||u||^/2|Au|i/2|u;|, for every u e K w e D(A), and w e H. (4.13) 
See e.g. [TemOl] for further details. 
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We next introduce the main functional spaces and operators associated with the memory terms in (4.9a), 
(4.9b). For any m > we define the Hilbert spaces 

Mm ■■= L^(K+,/i;i:»(A™/2)) ^ . is (R+,i:»(A"/2)) measurable,^ fi{s)\r]{s)\l^ds < ooj , (4.14) 

which are equipped with the inner products, 

[77,p]„:= / ^,{s){'n,p)mds= / ti{s) / • V dx ds (4.15) 

Jo Jo JV 

where A denotes the Stokes operator introduced above according to (4.10) along with its corresponding 
fractional powers recalled there. For notational convenience and in accordance with the classical notations 
for the Navier-Stokes equations we shall frequently denote [•,•] = [•,•]! and [[•,•]] = [•,-]2, which are the 
appropriate norms for 'weak' and 'strong' estimates respectively. See Remark 4.3 below. 
In view of term involving —ds appearing in (4.9b) we next take, again for any m > 0, 

Mm ■■= i?o'(K+,Ai;^(^"/')) = {v^Mm: d,ri G Mm^viO) = 0} . (4.16) 

On Aim we define an (unbounded) operator T according to 

Tt] = -9,77, 77 e M„(= Dom(T)). (4.17) 

We now summarize some basic properties of T used in the sequel. See e.g. [GP02, CPS06, CGHb] for further 
details. 

Lemma 4.1. Assume that the memory kernel /i G L^(M+) is nonnegative, non-increasing and satisfies (4.6). 
Fix 771 > and define T according to (4.17). Then: 

(i) For any for ij G Afm- 

[T77, 77]™<0. (4.18) 

(a) Suppose that rjQ G Mm o,nd that ^ G ^^^^^([0, cxj); Z?(A'"/^)). Then there exists a unique mild solution 
rj G C {[Q,oo); Mm) of the system: 

dtv' = Tr,*+at), V°=Vo- (4.19) 
Moreover rj has explicit representation (cf. (4. 8) J; 

[Vo[s - t) + £,[a)da if s > t, 

which is valid for any t > 0. 

Remark 4.2. In the case where K ~ I in (4.6) it is not hard to show that [Trj, rj] < —6/2. This significantly 
simplifies the estimates below in e.g. (4.26), (4.36). Indeed, it is only with the functional described in 
Lemma ^.2 (see [GMPZ08]) that the needed dissipativity properties for the memory variable rj can be achieved. 

We next define the product spaces which will serve as the phase spaces associated with the system (4.9). 
For m > let Jim ■= D{A"^/^) x Mm+i which we endow with the product norm: |j {u, 77)||^„ := + M^n+j. 
We will sometimes write x = (u, rj) G T-Lm and denote the underlying projection operators by Px ~ u and 
Qx — rj. In view of the above conventions, to simplify notations we will we will use often % := T-Lq and 
V := Ui and let 

\\{u,^)\\l = \u\' + W^ and 11(7., 77)11? = Wuf + Mf. (4.21) 
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Remark 4.3. The additional degree of 'regularity' in the memory space is dictated by basic the structure of 
system (4.9). See, for example, the estimates (4.25), (4.39) in the proof of Proposition 4-1 below. 

For the each of the spaces Mm we do not have the compact embedding of Mm+i ^i^to Mm- Thus, in 
contrast to H and V , V <$%. Of course this introduces additional complications for the proof of the existence 
of a global attractor in the extended phase for (4.9) and leads to the introduction of still further spaces. See 
Lemma 4. 3 below. 



With these definitions in hand we may now eapture the memory term in (4.9) a bounded operator 
M :Mra^ D(A"'/2) defined, for m > via, 



This definition is justified since, of. (4.6), fjL{s)ri{s)ds\^ < fj,{s)^/'^ n{s)^/'^ \ri{s)\,.ds < C^[»7]fe. 

With aU of these preliminaries in hand, using (4.10), (4.11), (4.17), (4.19) we now rewrite the system 
(4.9) in the abstract form 



where F = Pnf- As in [GGP05] we may now associate a dynamical system {S{t)}t>a with (4.23) as follows: 

Proposition 4.1. Suppose that F Cz V' and xq ~ {uo,r]Q) £ H. Then there exists a unique x — {u,r]) G 
C([0, 00),"^) with x{0) = Xq, such that Px = u E Lfg^{[0,oo),V) and satisfying (4.23) for every t > (in 
the appropriate weak sense). By taking, for any t > 0, S(t)xQ = x{t) = {u{t),rf) we may thus define a 
continuous semigroup {5(t)}f>o on the phase space X = % (endowed with its usual topology). 

This Proposition may be established along classical lines with a Galerkin approximation. See e.g. [TemOl] 
and additionally [GGP99, GGP04] for further details on Galerkin schemes in the history space context. 

4.1.3 The global attractor, invariant measures and related estimates 

With this basic framework now in place we proceed to prove that the dynamical system {5(t)}t>o defined 
by (4.23), in the sense made precise by Proposition 4.1, possesses a global attractor so that we may with no 
further efforts apply the results in Section 2 to (4.23). More precisely we prove: 

Theorem 4.1. Suppose in addition to the conditions imposed in Proposition 4-1 that F E H . Then the 
dynamical system {S{t)}t>o defined on % by (4.23) according to Proposition 4-1 has a connected global 
attractor A and hence the results in Theorems 2.1, 2.2 hold for {S{t)}t>o with {X,d) = {%, \\ ■ ||). 

The rest of this section is devoted to the proof of Theorem 4.1. The proof consists in verifying each of 
the two conditions imposed by Proposition 2.1. With the sufficient condition supplied by this Proposition in 
mind, we proceed in two step. At the first step we use energy type estimates to establish the existence of an 
absorbing set *8. Note that, in view of (4.18), the needed dissipation in variable 77 does not follow directly 
from these standard estimates. As such we employ an additional functional as dictated by Lemma 4.2 
introduced below. For the second step we achieve the splitting suggested by Proposition 2.1 by considering 
solution operators Si associated with no external forcing F and S2 with zero initial data. See (4.34), (4.35) 
respectively. Since, as mentioned in Remark 4.3, V is not compactly embedded in "H, we are required to 
introduce still further spaces at achieve the compactness desired for 5*2 . This is the significance of the space £ 
in (4.42) and Lemma 4.3 immediately following. Note that the formal estimates that follow may be rigorously 
justified in the context of a suitable Galerkin approximation scheme for (4.23). 




m • 



(4.22) 



dtu + vAu + M {Arj) + B{u) = F 

dtrj ~ Trj + u 

(«(0),7,0) = (7/o,%) 



(4.23) 
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Step 1: Energy estimates and dissipativity. By multiplying the first equation in (4.23) with u and 
making use use of the cancellation properties of B we infer 

i^|w|2 + i.||«||2 + [77,u]-(F,w). (4.24) 
We next multiply the second equation in (4.23) by Arj and deduce: 

~[v?^[Tv,ri] + [v,u]. (4.25) 
Adding (4.24) and (4.25), using (4.18) and making an obvious application of Young's inequality we have 

|(H2 + M2) + Hl^lP<a|F|^. (4.26) 

This inequality docs not evidence any dissipation in the 'memory variable' 77. We are therefore unable 
to directly invoke the Poincare inequality and the Gronwall lemma in order to conclude the existence of 
a bounded absorbing set as we would for e.g. the classical 2-D Navier-Stokes equations. See however 
Remark 4.2 above. 

To overcome this difficulty we rccaU the following Lemma (see [GMPZ08] and also [CDPGHPll, CGHb]) 
which we will use at this step and later on in Step 2 when we introduce a splitting of S(t) (cf. Proposition 
2.1, above). 

Lemma 4.2. Suppose that fi and k satisfy the conditions imposed in (4.5), (4.6). Fix any m > and assume 
that 770 G Mm, £. G LjociiO: 00); D{A"'/'^)) and that rj is the unique solution of (4.19) corresponding to ^ and 
r]o as guaranteed by Lemma 4-1, (H)- Define, the functional 

/•oo 

r™(0 - r,„(7y*) := / K{s)\rj\s)\',^ds. (4.27) 
Jo 

Then T„i (t) satisfies the following differential inequality, 

+ ^(r™W + PWfm) < 2/32|1mI|lhm+)ICWI™, (4-28) 

where j3 is the constant arising in (4.7). 

We now introduce, for each A > 0, the following functional, 

$^(t):=ri(77*) + A(Kt)P + [,7*]2). (4.29) 
Applying Lemma 4.1 with fc = 1 to the second equation in (4.23) we find that 

where :— 2/3^||/i||£,i(-g+), as arising in (4.28). Combining this observation with (4.26), we estimate 

1*° + ^(ri + m') < iC, - Ai^)\\ur + AC.\F\l,. (4.30) 

■ AXiiy V AXiu V ' ^' 4max{;3,A} ^ ' 



Take 
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where we recall that Ai the constant in the Poincare inequality and /3 comes from (4.7). With this choice of 
A and associated 7, we observe that 

1 



4/3 

and so with (4.30) and (4.31) 



With the semigroup notation x{t) = S(t)xo, xq = (uq, tjq) we therefore estimate using the Gronwall inequality 
that 



(4.32) 



Note that we have used (4.7) with (4.29) to infer the last inequality. We may hence infer the existence of 
an absorbing ball *Bo in H whose radius depend only on \F\y' and the universal C = Ci,^\-^^f^ in the final 
inequality above. More explicitly, (4.32) shows that, given any K > 0, there exists a time = ^1(^5^) such 
that for t > tl: 

S{t){xo e n : \\xo\\ < K} C So- (4.33) 

This is the first item required for Proposition 2.1. We turn to the second step. 

Step 2: The splitting property. We next exhibit a splitting of S{t) = Si{t) + 5*2 (t) taking the form 
required by Proposition 2.1. Given arbitrary initial data xq — (wqi^^o) G H we define {ui{t),'i]\) :~ Si{t)xQ 
to be the solution at time t of: 

dtui + vAui + M [A-qi) = —B{u, ui), 

dtm^Trji+ui, (4.34) 
(ui(0),r/J) = {uo,Vo)- 

We take {u2{t),r][) ~ S2{t)xo to be the solution of: 

dtU2 + vAu2 + M (Ar/2) ^ F - B{u, U2), 

dtm = Tin + U2, (4.35) 
(7/2(0),%°) -(0,0). 

Note that the u appearing in the nonlinear terms of both (4.34), (4.35) is the first component of the solution 
of (4.23); in other words u = PS{t)xo- Clearly, given any xq £ ^, we have S{t)xo = Si{t)xo + S2{t)xo- 

The estimates for 5*1, in view of Proposition 2.1, (ii) are carried in Ji and are very similar to those in 
Step 1 for S. Indeed following the arguments leading up to (4.26) we obtain 

|(KP + NiF) + Hl^ilP<0. (4.36) 

We now make a second apphcation of Lemma 4.2. Defining, for A > 0, ^\{t) :— ri(?y^) + A(|wi(t)p + [77^]^), 
arguing as above in (4.30) and tuning A and 7 as in (4.31) we infer 

\\^iml<C.MAe~''\M\l), (4.37) 
where xi{t) — Si{t)xo. We infer that for every K > 0: 

sup \\Si{t)xo\\o 0. (4.38) 

\M\o<K 

so that indeed 5*1 plays the desired role required for Proposition 2.1, (ii). 

In remains to establish the requirement on S2 given in Proposition 2.1. (ii); we must show that, for each 
t > 0, and each bound subset B of H that S2{t)B is a precompact subset of "H. To this end we first show 
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that S{t)B is a bounded subset of V and then, owing to the coniphcation that V is not compactly embedded 
in H, we must take further steps. See Remark 4.3 and Lemma 4.3 below. 

In order to have an equation for 1 1— |j5'2(t)a;o||i we multiply the first equation of (4.35) by Au2 and the 
second equation by A^r]2- With (4.18) one finds that 

|(!|M2f + M') + iy\Au2\' < 2\{B{u,U2),Au2)\ + C\F\\ 
By applying standard estimates on the nonlinear terms, (4.13), we find: 

2\{Biu,U2).Au2)\<C\u\'^'\\u\\'^^U2\\'/Mu2\'^' <CM^h^^^ 

Hence: 

|(lk2f + [[V2]f) < C.i\uf\\uf\\u2f + \Ff) < C4\u\'\\uf{\\u2f + [fef) + |F|^). (4.39) 

Let T(t) ~ Ci' /o [""PlI^lP where Ci, is precisely the constant appearing in (4.39). Returning to (4.26) we 
may observe that for any xq = {uo,r]o), we have supj,<( |up + i' \\u\\ds < |uoP + tC,j\F\v', so that clearly 
e"*"^*'' < oo for each t > 0. Applying this functional to (4.39), integrating in time appropriately, and noting 
that S'2(0)a;o = we find that 

sup \\S2{t)xo\\i = sup {\\u2ir)f + M]?) < f e'^'^-^^'^\F\^ds < C,,t\F\^. (4.40) 

0<r<t 0<r<t Jo 

This shows that, for every bounded set B C H and for every t > that S{t)B is a bounded subset of 
V. However, as described above in Remark 4.3, V is not compactly embedded in H and further steps are 
therefore required to established the desired compactness of 52. 

To compensate for this difficulty we follow previous works (see e.g. [GGP05, CPS06]) and introduce some 
additional spaces. On A^i we define 'tail functional' 

T^(77) =sup a - / ^^{s)Us)fds, (4.41) 

'^>l J(0,1/(t)U(<t,oo) 

and consider the subspace 

£ = {veAf2:T^ir]) <oo} cAf2, (4.42) 

which we endow with the norm, |?/|| :~ [[?/]]^ + [T'q]'^ + T^(r/). Under these definitions £ may be shown to 
be a closed subset of A^2 (relative to | • Is) and hence is a Banach space; see [PZOl]. We thus define the 
product space V = V x £ G V and endow V with the norm, 

\\{u,v)\\l ■■= M' + Hl = \\{u,vm + [Tv? + 1,{v)- (4.43) 

We have the following compactness results and related estimates on the evolution equation for the history 
variable. See e.g. [GGP05, CPS06, CGHb]. 

Lemma 4.3. Assume that the memory kernel /i G L^{M.'^) is nonnegative, non-increasing and satisfies (4.6). 
(i) £ is compactly embedded in Aii and hence V is compactly embedded in %. 

(ii) Suppose that rjQ ~ {f' , and that £ {{0,oo);V). Letrf be the unique mild solution in C{[0,oo),Aii) 
of, 

dtrf ^Ttj' +^{t), 7/0 = 0, (4.44) 

^With some modifications we may assume merely that rio a £. Since this is unncedcd for our current purposes we do not 
state the Lemma in this greater generaUty here. 
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guaranteed by Lemma 4-1, (in). Then there exists a universal positive constant C ~ C'^ (depending on 
fj. but independent oft and the data) such that 

[r,7f + 1^{7f) < C sup UisW < oo. (4.45) 
se[o,t] 

With Lemma 4.3 in hand and having aheady estabhshed (4.40) we finally proceed to show that S2 is 
compact. Fix any xq = (uo,77o) G "H. Observe that, due to (4.40) we may infer that 

sup \\PS2{t)xof < sup \\S2{t)xo\\l := J{t) < 00 (4.46) 

0<r<t 0<r<t 

where we note that J{t) is independent of xq. We now apply Lemma 4.3, (ii) with ^ = PS2{t)xo in (4.44). 
From (4.35), (4.45) and (4.46) we infer 

[TQS2it)xo]'' + 1^,iQS2it)xo) < CJ{t) (4.47) 

where C is the constant appearing in (4.45) and is independent of xq. Combining (4.46), (4.47) and cf. 
(4.43) 

sup \\S2{t)x4l<CJ{t) <<x, 

\\M\a<K 

for any K > 0. In other words we infer that for every bounded set B C Ti and for every t >Q that 6*2 (t)i? is a 
bounded subset of V C V and so conclude from Lemma 4.3, (i) that 6*2 in compact. The second requirement 
for Proposition 2.1 has therefore been fulfilled and the proof of Proposition 4.1 is finally complete. 

4.2 Neutral delay differential equations 

A retarded functional differential equation (RFDE) describes a system where the rate of change of its state is 
determined by the present and the past states of the system. If, additionally, the rate of change of the state 
depends on the rate of change of the state in the past, the system is called a neutral functional differential 
equation (NFDE). When only discrete values of the past have infiuence on the present rate of change of 
the system's state, the corresponding mathematical model is either a delay differential equation (DDE) or a 
neutral differential equation (NDDE) . The theory of RFDEs and NFDEs is both of theoretical and practical 
interest, as these types of models provide a powerful framework used in the study of many phenomena in 
the applied sciences; for example in physics, biology, economics, and control theory to name a few. 

The initial development of neutral functional differential equations has been intimately related to some 
particular transmission line problems modeled by hyperbolic PDEs such as the telegrapher's equations. In 
the seminal works [AM60], [BM64], it was observed that certain linear hyperbolic PDEs with (non)linear 
boundary conditions were equivalent to an equation of the form xt — 2^Xt = f{xt) where £ is an operator which 
does not depend upon the values of x{t) but only on values in the past. In other words the evolution equation 
of the system may be characterized by the simultaneous presence of delayed and non-delayed derivatives. 
We will explore this relationship in further detail for the specific case of the telegrapher's equation, with 
simplified, linear, boundary conditions (4.55), below in subsection 4.2.1. 

The linear boundary conditions (4.55) have been chosen here merely to give a simple illustration of how 
NDDEs may arise in problems modeled by hyperbolic PDEs. More complex, physically realistic boundary 
conditions, such as nonlinear or dynamic boundary conditions, lead to many other (nonlinear) NDDEs. See 
e.g. [IIL93, BC04, BKW07]. In [BC04] for example, a nonlinear NDDE model derived from the telegrapher's 
equation, but with nonlinear boundary conditions, was employed to model an actual electronic device. Here 
the delay may be seen to arise partially from the signal passage time through the transmission line itself. 
The presence of rich, chaotic delay-induced dynamics, was observed experimentally and shown to be in good 
agreement with the numerical simulation of their nonlinear NDDE model. More precisely it was shown 
that, in the appropriate regime, the attractor determined by this NDDE presented a similar coarse-grained 
structure in comparison to the attractor reconstructed from the empirical observation of the electronic device 
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they were studying. This particular NDDE has been further analyzed in [BKW07]. Here it was shown that 
homoclinic bifurcations were at the origin of the complicated dynamics observed in [BC04]. 

Of course, to determine if such complicated dynamics may be associated with the existence of a "com- 
plicated" invariant probability measure is out of reach due to the infinite dimensional character of the phase 
space X = C([— r, 0], M").^ Indeed, in view of the rich chaotic dynamics observed in [BC04, BKW07], even 
the question of the existence of an invariant probability measure for such systems seems to be nontrivial 
and, to the best of the authors' knowledge, this question seems to have never been previously addressed in 
the literature. 

Theorem 4.2 below yields a general, abstract result concerning the existence of invariant probability 
measures of NFDEs and furthermore associates these measures with certain temporal averages. For each 
initial Borel probability measure mo, we obtain an invariant measure m as a generalized Banach limit of mo 
through the time average of the semigroup generated by the NDDE. Note that, by our construction, the 
support of any invariant measures obtained by Theorem 4.2 (see also Theorem 4.3) is furthermore contained 
in A. This does not preclude the possibility that the invariant measure m could be supported by the whole 
of the global attractor with its complex geometry. Our results apply to any NDDE possessing a global 
attractor A; cf. Theorem 2.1, 2.2 and Definition 2.1 above. Note that, in this case, such a semigroup has 
no smoothing effects in finite time (i.e. the semigroup does not have a compact absorbing set) and acts on 
X = C([— r, 0]; M") which is a non-reflexive Banach space. As such, the results in [LRRll] do not apply 
for such equations. To provide concrete situations where NDDEs possess a global attractor. Theorem 4.3 
provides a novel sufficient and "checkable" condition which yields a broad class of interesting examples. 

The rest of this section is organized as follows. We begin by further detailing the relationship between 
NDDEs and hyperbolic PDEs with an example involving the telegrapher's equation (with simplified boundary 
conditions). In subsection 4.2.2 we briefly recall some preliminaries from the general theory of NDDEs which 
allow us to establish the main abstract result. Theorem 4.2, which establishes the existence of invariant 
probability measures for a broad class of NFDEs. In subsection 4.2.4 we formulate and prove a result of a 
more practical nature. Theorem 4.3, which gives concrete sufficient conditions for the results in Theorem 4.2. 
The final subsection applies Theorem 4.3 to some examples related to the classical NDDE model proposed 
in [BM64]. 



4.2.1 NDDEs and hyperbolic PDEs 

We can get some idea of how NDDEs (and also difference equations with retarded arguments may) be derived 
from certain classes of hyperbolic PDEs by considering an example involving the telegrapher's equation with 
particular linear boundary conditions in some details. We shall see that the type of retarded system derived 
depends on the boundary conditions. For this purpose, we recall the classical problem considered in [BM64] 
of the determination of the current / and voltage ^ in a transmission line terminated at each end by linear 
resistors or nonlinear circuit elements such as diodes. Note that we adopt here a different approach than 
the one classically presented in [HL93]. The reasons of this, which are related to stability and regularity 
issues, will be clarified in the course of the forthcoming discussions. In particular, the exposition adopted 
here allows us to illustrate an important relationship between difference operators and NDDEs. As we shall 
see, cf. Lemma 4.4 below, the structure of a certain difference operator associated with the type of NDDEs 
(or NFDEs) we consider here determines the non-compact nature of the associated semigroup. 

Let us take the x— axis in the direction of the line, with ends at a; = and x = 1. It is assumed that the 
voltage V and current / at each point along the line are governed by the following telegrapher's equations 

d^V ^ - Ldtl, 

4.48 

dj ^ - cdtV, ^ ' 



°This question seems intractable since there is no equivalent notion of the Lebesgue measure to employ as a universal 
reference measure mo. We thus have no means to declare that a measure m to be "complicated" in terms of a lack of absolute 
continuity with respect to a reference measure. 
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where L and C are respectively the inductance and the capacitance, both of which will be assumed to be 
constant. It is well-known that the solutions of (4.48) can be represented as a superposition of traveling 
waves moving towards x — {) and x = \ with velocity c = (LC)~^/^. We have then for x € (0, 1) and < > 0, 



V{x,t) =-[(j>{t - c~^x) + yj{t + c~^x)], 

H^, t) ~m - c-'x) - i,{t + c-'x)] = 

ZLc I 



(4.49) 



c ^x) ~ip{t + c ^x)], 



where (f> and ip are some unknown functions to be determined from the boundary conditions. See e.g. [Smi64, 
Section 164]. Frequently, (j) (resp. ip) is called a progressive wave (resp. reflected wave) since its argument 
is a solution of the increasing (resp. decreasing) characteristic dt/dx ~ c^^ (resp. dt/dx ~ — c^^); see e.g. 
[ZT86]. 

Let us consider first the following linear boundary conditions 

y(0, t) + i?o/(0, t) = E, I{l,t) = 0, for t > 0, (4.50) 

where i? is a voltage which is assumed to be constant for simplicity. Physically, this corresponds to a trans- 
mission line that is "disconnected" from the nonlinear circuit element (diode) at the right hand boundary, 
i.e. /(l,i) = 0, with some resistor of resistance Rq at the left boundary. Such boundary conditions are 
expressed in terms of 4> and '0 as 



(1 



il-RoVC/L)m=^E, 
: (f>{t-VLC). 



(4.51) 



Consider the Cauchy problem associated with (4.48) - (4.50) and given initial data V{0,-) = Vo(-) and 
/(O, •) = /o(')- From the theory of first order PDEs, we know that this Cauchy problem is well-posed provided 
that (x, t) belongs to the triangular domain T determined by the characteristics according to. 



1 



{{x,t) e [0,1] 



1 



X X 

- + - <t <- and 

c c c 



X " 



X 1 

- <t< 

c c c 



Now, in order to determine (j) and from the initial data and boundary conditions, we observe that since 
X £ [0, 1], can be determined from Vq and Iq only on [— a/LC*, 0] according to </'(— f ) ~ Vq{x) + L/CIq{x). 
Similarly we initially determine ijj only on [0, -\/ LC] via ipi"^) = Vo{x)— ^/lJCIo{x). It is then equation (4.51) 
which allows one to determine (j) and ip for all of R. To be more precise, let us introduce ip{t) = -tpit + y/LC); 
now since ip is known on [0, \/LC], "ip is determined on [— -s/LC, 0], as (j> is- We can therefore rewrite the 
system (4.51) as the following difference equation in continuous time: 







r 2_E 1 


+ 


1 + Rq^C/L 


i,{t - r)_ 








(4.52) 



where r = V LC. This delay r may therefore be interpreted as the characteristic time associated with wave 
propagation; it represents the transmission time through the transmission line itself. 
The system (4.52) may be written in a compact form as 



x{t) ~ Bx{t -T)^f, 



(4.53) 



where x{t) = [(j){t) ip(t)]'^ , B is the matrix arising in (4.52) and / denotes the constant vector in (4.52). By 
introducing the operator Dq^ := $(0) — B^{—t) acting on functions $ G X := C([— r, 0], R^) this equation 
takes the functional representation 

Doxt = /, (4.54) 
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where Xt G X, for all t and we take Xt{0) = x{t + 0) for all 9 E [— r, 0]. When supplemented with the 
initial data xq = ip]'^ , where (t>{0) = Vo{-c0) + yJTJChi-cQ) and ^{B) = Vq{cQ + 1) - s/LjCIo{ce + 1), 
with 9 g [— r, 0], equation (4.54) (along with (4.49)) gives thus another representation of the Cauchy problem 
associated with the original boundary value problem (BVP) (4.48)-(4.50). 

Remark 4.4. Let xq : [— r, 0] — > be some initial condition in X and define x(t), starting from xq according 
to the so-called method of steps (see e.g. [SmillJ); having determined x(t) on the interval [— t, (fc — 1)t) 
for k > 1, we determine x{t) on [{k — 1)t, fcr) using (4.53). Then a;(0"'') — x{0~) = Bxq{—t) + / — xq{Q) 
with a;(0"'"^~) = lim(_^Q+/- x{t); unless xq satisfies (4.53) at t ~ 0, we obtain a discontinuity in x at t ~ 0, 
i.e. a:(0'*') — a;(0^) 7^ 0, that will propagate to each t = kr, fc G N. To avoid such a phenomenon we 
need to assume that Bxo{~t) + / = a;o(0) which translates to 1^(0,0) + i?o^(0,0) = Vo(0) + i?o^o(0) = E 
and 7(1,0) = /o(l) — in the original problem BVP (4.48)-(4.50), requiring in other words that the initial 
condition Iq, Vq must satisfy the boundary conditions (4.50). This phenomenon is reminiscent of the well- 
known phenomena of singularity propagation along the characteristics, typical of hyperbolic PDEs, see e.g. 
[ZT86]. 

Now, if we consider instead of (4.50) the following linear boundary conditions 

^^+Ro^=E, ^ = 0, (4.55) 
dt dt ^ dt ' 

similar computations lead us to the following functional equation which is now an NDDE of the form 

j^D^x^ = /, (4.56) 

giving another formulation of the Cauchy problem associated with the boundary value problem (BVP) 
(4.48)-(4.55). Note that, here the operator Dq is still defined as Dq^ ^{^) ~ B^{—t) acting on functions 
$ S X := C([— r, 0], M^). From this introductory example, we can sec how the difference operator Dq arises 
naturally in the NDDE formulation. Further details concerning the relationship between difference equations 
and NDDEs, will be given in the next section. 

Remark 4.5. As mentioned above, the telegrapher's equations are not the only examples of hyperbolic PDEs 
which can be reformulated as systems involving retarded arguments of similar nature to NDDEs. For instance, 
in modeling of the El Nino-Southern Oscillation (ENSO), the authors in [GTOOJ have proposed such a neutral 
formulation derived from a model introduced in [JN93a, JN93b]. In these later works the ENSO dynamics 
were described by a linearized shallow-water (hyperbolic) model on an equatorial (3-plane coupled with an 
advection equation describing the sea surface temperature changes at the earth's equator. In [GTOO], by 
dropping the spatial dependence in the advection equation, the authors reformulated the resulting PDE-ODE 
model as a system made up of difference equations coupled with ODEs. Note that, coupled systems of ODEs 
with difference equations arise in other physical applications dealing with wave propagations phenomena, and 
are strongly related to the theory of NDDEs, see [NicOl ]. 



4.2.2 Existence of invariant measures for NFDEs: preliminaries and abstract result 

We recall here some preliminaries concerning neutral functional differential equations used in the statement 
and the proof of Theorem 4.2 below concerning the existence of global attractors and hence of invariant 
measures for the semigroup generated by such equations. In the next subsection (see Theorem 4.3 below) 
we will provide concrete conditions which allow us to apply this abstract result to a broad class of NDDEs. 
While, as we said above, the existence of invariant measures for NFDEs is new, the exposition here largely 
follows the framework presented in [HL93] , to which we refer the reader for further details on the extensive 
general theory. 

Throughout what follows we take (•, •) be the inner product on R", with | • | the corresponding norm. Let 
T > 0, and take X := C([— r, 0],R") be the space of continuous functions taking [— r, 0] into M" which will 
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be endowed in what follows with the topology induced by the supremum norm viz. 

\4>U ■■= sup 10(0)1, for ex. 

-T<e<Q 

We denote by A^„(K), the space of ?i x n matrix with real coefficients which we endow with the natural 
norm i.e. := sup{|Afx| : x G M", \x\ = f} for any M e A1„(M). 

Throughout this section we consider a fixed map M taking [— r, 0] into A^„(M) which is of bounded 
variations in the sense that, 

Var(A/) sup^ ~ M{e,+i)\\ < oo, 

l-r.o] 

the supremum being taken over all finite partitions — r ~ Oq < 9i < ■ ■ ■ < 9^ = oi the interval [— r, 0]. We 
assume furthermore that 

Vse (0,t], Var(Af) < 7(5), 

[-s,0] 

with 7 a continuous nondecrcasing scalar function for s > 0, such that 7(0) ~ 0. In language of measure 
theory, this last requirement implies that the vector measure M (after identifying M„ with R" ) is nonatomic 
at zero, i.e. M(0) — M(0~) = 0. In other words AI does not attribute a non-zero measure to {0}. This 
assumption about M being non- atomic at zero gives a sufficient condition to develop a theory of existence 
and uniqueness of solutions of NFDEs of the form (4.58); see [HL93] and Proposition 4.2 below. 
With this M in hand we define a bounded linear operator D : X ^ M" according to 

r" 

Dcj) 0(0) - j [dM{e)](f){9). 7 (4.57) 

In (4.57), the notation dM{9) before the integrand (j){9) emphasizes that M{9) is a matrix and 4){9) is a 
column vector so that the integral is column vector- valued. 

Remark 4.6. A classical and fundamental variety 0/ examples 0/ such operators D arises when we take M 
to he a step function. This leads to 

N 

i?0 = 0(O)-^Bfc</)(-Tfc), 
fe=i 

where the Bk 's are n X n matrices and the delays satisfy < ti < ... < r^r < r. This type of operator D is 
very often encountered in the applications, see e.g. [MN07] and the references therein. 

Let X be an R" valued, continuous function defined on an interval [— r, T] for some T > 0. For t E [0, T], 
we define Xt G X as being the "copy" of x over the time interval [t — T,t] shifted down to [— r, 0], i.e. for 
each t e [0,T], we define Xt & X via Xt{9) = x{t + 9) for all 9 e [— t, 0]. For simplicity, let us assume that 
/ : X — > R" is as regular as is needed for a moment. A neutral functional differential equation (NFDE) is 
then given by the following relation: 

j^Dxt = f[xt), (4.58) 

where ^ denotes the right hand side derivative at t. The initial data is given as element in e X. To make 
sense to (4.58) it is required that t — > Dxt is as a map with values in R"; cf. [HL93]. More precisely we 
have the following definition of solutions of (4.58): 

'^Note that, by application of the Bartle-Dunford-Schwartz theorem (see [BDS55, Theorem 3.2]), any bounded linear operator 
L from X to R" may be represented as Stieltjes integral f^.j.[dN{9)]4i(d) where N is an valued measure on Borel sets of 
[— T, 0] which is of bounded variations. Operators having the form of (4.57) may therefore be interpreted as those which have 
the representation corresponding to A' = 5o ~ ^ with M being non-atomic at zero and <5o being the Dirac measure at 0. 
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Definition 4.1. For a given (j) € X, we say that Xt{-',4>) is a solution of (4.58) on the interval [0,a^), 
> 0, with initial data (j) at t = 0, if xo{-, (jy) = 4>{-), Xt{-', 4>) G X , for all t G [0, a^), the map t — > Dxt{-] 4>) 
is in C^([0, a0),]R"), and Xt{-,(t>) satisfies (4.58) for all t e [0, a^). 

Remark 4.7. 

(i) Note that, since M is assumed to he non atomic at zero, we have that J_^[dM{e)]<j){9) for 

any (j) €z X . In this limiting case Dxt ~ x{t) for all t > 0, which means that the the case of ODEs 
corresponds formally to t = with the well-posed problem of x ~ f{x) reducing from Definition 4-i to 
the classical sense of Hadamard. 

(a) Similarly to the case of difference equations, starting from initial data, singularities may propagate 
but at the level of the derivatives (cf. Remark 4-4 ) o,nd we do not have in general that t xt is hut 
in such a case, (4.58) reduces to Dxt = f{xt), an equation in which the derivative occur with delayed 
arguments; see [HL93] for further details. 

We recall now the main results on existence, uniqueness and continuous dependence of the solutions of 
(4.58); and wc refer to [HL93] for a proof of these properties. 

Proposition 4.2. Let D be as given by (4.57) and consider any f be in C^(X, R"). Then, for any (p £ X , 

there exists a unique solution Xt{-;(t>) of (4.58) through (f> at t = defined on a maximal interval [0,a^), 
such that if < oo then lim |xt(-;0)|oo = oo. Furthermore, the map {t,(j)) — > Xt{-](j)) is continuous from 

the set [0, a^) x X into X. If we assume that all solutions of (4.58) are global in time, i.e. we assume that 
— oo for all (j) £ X , and we introduce SDj{t)4> '■= xt{-,(j)), then {SDj(t)}t>o is a continuous semigroup^ 
acting on X. 

Hereafter in this section, we will assume that D as given by (4.57) and / € C^{X,M."-) are such that all 
solutions of (4.58) are global in time.^ 

We next describe a decomposition of the semigroup Snj'it), given in precise form in Lemma 4.4 below, 
which splits Snjit) into two semigroups, one of which is compact and the other of which decays exponentially 
towards zero. As we shall see, this decomposition will depend on the asymptotic stability of the zero 
solution of a certain difference equation associated with D; cf. (4.59), (4.60) below. Note that, in view of 
Proposition 2.1, Lemma 4.4 will provide a sufficient condition for the existence of a global attractor. 

To formulate this result we will assume that the operator D is given by, 

D<j> = Dq4) - j M{9)(f>{9)d9, with, 

(4.59) 

Docjj = (/.(O) - Bk(bi-Tk) for <peX, 

k=l 

where the Bk's are n x n matrices with the delays of the form < ri < • • • < rjv < ''"7 and where the 
family of matrices {A/(0)}Qg[_.r,o] satisfies \\M{9)\\d9 < 00. We now consider a linear difference equation 
associated with the operator Do given by: 

N 

Doyt = y{t) - J2 ^^y^^ ~ ^fc) = 0' ^ ^ 0. (4.60) 

Clearly, any initial value yo of this equation must live in the null space Af^Do), so wc can restrict our attention 
to Xug the closed subset of X defined according to: 

Xdo := ^ n MiDa) = {0 e X : Do^ ^ 0}. 

^In the sense recalled at the beginning of Section 2. 

^See Theorem 2.1 in [HC69] and references therein, for classical conditions ensuring existence of global solutions in time. 
See also Lemma 4.5 proved below. 
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Then, for any (p G ^D^n it can be shown that there is a unique solution t — > yt{-;(j)) of (4.60) which is 
continuous from [— r, oo) into X^q and coincides with (p on [— t, 0]; see [HL93].^° This imphes that the 
translation along the solutions of (4.60), 

defines a C" semigroup on X^^. As shown in [Hal77, Theorem 3.3. p. 284], the infinitesimal generator Aq 
of this semigroup is given by 

Ao(j) = 4> for ^ e V{Aa) = {0 € Xd,, : e ^Do}, 
and the spectrum of Aq, <7{Ao), is then given by 

a{Ao) = {A e C : det(H(A)) = 0} 

where, 



N 



k=l 

This last relationship can be heuristically understood by looking for nontrivial solutions of (4.60) of the form 
VtiG) = e^^*~^^^vx, 6 € [— r, 0] where vx is some vector living in ker(_ff (A)). 

As we already mentioned, Lemma 4.4 below is conditioned to the stability of the zero solution of (4.60). 
This stability is related to the location of the spectrum of Aq with respect to the imaginary axis, which is 
itself determined by the roots of the characteristic equation dct{H{X)) = 0; see e.g. [MN07]. According to 
[Hal77, Theorem 4.1, p. 287, Lemma 3.3, p. 284], we deduce that if the real part of the rightmost eigenvalue 
of ^0 is located in the left half-complex plane, i.e. if 



^y^Ao ■= sup < Re(A) : det 



N 



k=l 



= ^ < 0, (4.61) 



then is exponentially stable in Xoo for Sooit) in the sense that there exist C > such that for alH > 
and all cj) G Xb^ we have IS'^^ (i)0|oo < Ce^"*|(/)|oo, with a ~ '^Aa- Ii^ what follows we will summarize this 
property by just saying that Do is stable. 

With this background now in place, the desired decomposition lemma may now be stated as follows: 

Lemma 4.4. Consider D, and f such that (4.58) is globally well posed in the sense of Proposition J^.Z. 
Let {5'£)./(<)}t>o be the continuous semigroup generated by (4.58) with D satisfying (4.59) such that the 
difference operator Dq is stable in the sense that ^Kao O'S given by (4.61) is strictly less than zero. Then there 
exists a (time-independent) linear bounded operator : X ^Dq such that, 

SDj{t)^SDo{t)o^ + UD{t), t>0, (4.62) 

where {UD{i)}t>o is a compact semigroup on X, and for every K > 0, 

sup \SDo{t) ° '^i4')\oo *^°°> exponentially with uniform decay rate ^Raq- 

Proofs of this lemma may be found in [HC69, CK70]; cf. also [HL93] for a unified treatment. Note that 
the noncompact part of Sojit) comes from the semigroup Sooit) which is associated with the difference 
equation (4.60) for which = is solution. The linear bounded operator ^ arises essentially in order 
to map X into Xdo, the domain of Sooit)- Consider, for instance, the case M = 0, and Bk ~ for all 
k G {1, ...jN} in (4.59), i.e. when the NFDE (4.58) becomes a standard RFDE with no retarded arguments 



^''Notc that such a yt can be derived from 4> G Xjj using the method of steps; cf. Remark 4.4, above. 
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on the derivative. Here, ^' is simply defined as the shift ^((/)) ~ 4> — 0(0) so that '5(0) G Xjjg where, in this 
case, Xdo = {(p E X : 0(0) =0}. In the general case, Vl' = Idx - ^ o Dq where $ : K" Xdq is the right 
inverse of Dq, i.e. Dq o $ = Idjjn. Indeed, for such a VP, we get ^"(0) = — $ o Dq^c/)) for all S X which 
gives when we compose to the left by Dq so that, as desired, 4'(0) G X^ij, for all G X. The existence of 
such a right inverse $ is established in [Hal??] for operators D of the form (4.57) with M nonatomic and in 
particular for operators such as Dq given in (4.59). 

Combining Lemma 4.4 and Proposition 2.1, we conclude the following general result concerning the 
existence of a global attractor for a wide class of NFDEs. The existence of invariant probability measures, 
associated with temporal averages via generalized limits, then follow for such NFDEs in view of Theorem 
2.1 and Theorem 2.2. 

Theorem 4.2. Let X := C([— r, 0], K"). Assume that D satisfies (4.59) with Dq stable in the sense that 
'^Kaq < 0, where is defined according to (4.61). Suppose that for this D and some given f G C^{X,M") 
that all of the solutions of the associated system (4.58) are well-defined for all t > in the sense of Def- 
inition 4.1. Let {SD.f{t)}t>a be the continuous semigroup generated by (4.58) on X. If there exists an 
absorbing set *B C X, for So.f then the omega-limit set, w(*8), is the global attractor of Sdj- Moreover, 
for any given generalized Banach Limit LIM , and for any mo G Pr{X) there exists an invariant measure 

T->oo 

m G Pr{X) for {SDj{i)}t>o whose support is contained in A and such that 

[ ip{u)dm{u) = LIM- [ [ (p(5d /(t)u)dmo(M)dt, for any ip € Cb{X) . (4.63) 

Jx T^ooT Jo Jx 

Furthermore, if mo ~ 6^, the Dirac measures for some G X , then (4.63) holds for any if G C{X). 

Proof. The result follow as an obvious consequence of Lemma 4.4, Proposition 2.1, Theorem 2.1 and Theorem 
2.2. □ 



4.2.3 Invariant measures for NDDEs 

In view of its generality, Theorem 4.2 provides a powerful tool that may be employed to establish the existence 
of a global attractor and hence of invariant probability measures for particular NDDEs. 

Of course, in practice, we need to verify the condition UIao < 0, which ensures the splitting of the 
semigroup generated by the NDDE of interest, and furthermore to establish the existence of a bounded 
absorbing set. We next provide below a useful criterium for the verification of this latter dissipativity 
condition. This criteria is satisfied for a special but still rather broad class of NDDEs. The verification of 
9\ao < has been and is still an intensive topic of research. Many criterion exist depending on the situation 
of interest; see e.g. [MN07]. In what follows we will need only the Schur-Cohn criterion, see e.g. [NicOl]. 

Lemma 4.5. Letr be apositive constant, and B anxn matrix with real entries. Let g be in C"'^(R" xM"; M"). 
Consider the neutral delay differential equation given by 

^ (xit) - Bxit - r)) = .g(.T(<), a:(i - r)), for t > 0, (4.64) 

on the phase space X := C([— t, 0], M"). 

Assume that there exist a,7 > 0, /3 G K. such that 

{u - Bv, g{u, u)) < 7 - a|up + I3\v\^, V u, w G M", (4.65) 

and that. 



/ 1 — e~"'^ 

£:= ||B|| + \/(l + ||B||)2e-"- + 2(/3 + a||B||2) < 1. (4.66) 

V a 

Then the solutions of (4.64) are global in time, and the ball B{0, 2r-\/^^Q C'^) in X ^ C([—t, 0],R") with 
^ ~ V '^''^^^a absorbing for the continuous semigroup generated by (4.64) on C([— r, 0], M"). 
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Proof. Equation (4.64) may be written in the functional form -^DoXt = f{xt) with I3o0 '■= 0(0) — 
and /(0) := 5(0(0), 0(-t)) for any (/) G X, and since g G C^(R" x R"; R"), we have that / € C^{X, M"). The 
local existence and uniqueness of solutions up to a maximal time (4.64) is ensured by Proposition 4.2. 
From the same proposition, in order to conclude to the global existence of solutions of (4.64), i.e. to show that 
Q!0 = oo for every (j) € X, it suffices to show that for any </) and for any T > 0, supjgjo,mm{T,o!^}) I^t(s 'P)\oo < 
oo. The estimates provided below for the existence of an absorbing ball contain such an argument implicitly 
and existence of a continuous semigroup associated to (4.64) follows from Proposition 4.2. We turn next to 
these estimates. 

We first show that, for any < < ^2, 

\x{t2) - Bx{t2 - t)|2 <e-"(*=-*i)|.T(ti) - Bx{h - t)P + 2^(1 - e-"(*^-*i))+ 

2(/3 + a||Bf) sup |a;(s)p(l-e-"(*=-*i))/a. ^'^'^'^^ 

sS[ti— T,t2— r] 

This inequality may be derived by considering the quantity u(t) := |a;(t) — Bx{t — t)^ = {x(t) — Bx{t — 
T),x{t) - Bx{t - t)). Since Dqx G C1([0, c»), R"), u is diffcrentiable we find with (4.66) that 

«'(i) <2(7-a|x(i)P+/3|x(t-T)|2), 

Now since 2|x(i)p > u{t) - 2\Bx{t - t)P, we infer 

u'{t) < 27 - au{t) + 2(^ + a\\Bf)\x{t - t)^, for t > 0. 

The inequality (4.67) is then easily derived by multiplying this last inequality by e"* and integrating between 
ti and t2- 

Let us take now t to be in (0, r]. We use (4.67) with ii = and t2 = t, it follows that, 

\x{t) Bx{t r)p < e-"*|x(0) - Bx{-t)\^ + 2^(1 - e""*) + 2(/3 + a|15f)|0|2 li^^, 

a a 
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and since |a;(0) - Bx{-t)Y < (1 + \\B\\Y\(t)\l^, we get 



\x{t)~Bx{t-T)\' < (l + \\B\\f + 2{P + a\\Bf)^^^^ |0|^^ + 2^(1 - e-"-) 



a 

which gives finally. 



a 



|a;(i)| < Col^loo+r, forte(0,T], (4.68) 



with £0 := \/(l + ||B||)2 + 2(/3 + a||B||2)ii:i^-^ + and r = ^2^(1 - e-"^). For i > r we take 

ti ^ t — T and t2 = t in (4.67). By arguing in similar manner to the previous case we obtain: 

\x{t) - Bx{t - r)|2 < [(1 + llSlD^e-"- + 2(/3 + a\\Bf)^-^—^^] \xt-r\lo + r\ (4.69) 
which leads to, 

\x{t)\ < €\xt-^r\oo + r, for t > t, (4.70) 
with £ given in (4.66). We infer from (4.68), (4.69) and a simple induction that, for every fc > 0, 

fe 

|a;(i)l < e:*^Co|</>|oo +r-^e:^ forte (fcr,(fc + l)T]. (4.71) 

The existence of the bounded absorbing set for the semigroup generated by (4.64) therefore follows, com- 
pleting the proof. □ 
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Note that the dissipation condition (4.66) is vaHd for "large" delays provided that 2/3 < a and is 
appropriately chosen with respect to ^. More precisely we have the following: 



Lemma 4.6. (Dissipation for large delays) Suppose that 2(3 < a and \\B\\ < y2{l ^ §) ^ 1- Then (4.66) 
holds for all T > T*, with t* = -i log ( ■ 

Proof. Direct manipulations show that (4.66) is equivalent to: 

P(i|i3||)e-"^<P(|li3||+2), (4.72) 

where P is the polynomial function given by P{x) = —[(a; — 1)^ + 2(^ — 1)]. Under the assumed conditions 
on /3 and a, 




and P{x + 2) possesses two real distinct roots, namely = — 1 — \/2(l — ^) < 0, a;2 = — 1 + y 2(1 — ^) > 
and reaches its maximum at —1. As such, for all x G 



0, -y/2(l - §) - 1 j , it is clear that P{x),P{x + 2) > 

and also that P{x) — P{x + 2) is strictly increasing and non-negative. Since, by assumption, falls in this 
range, (4.72) holds for all r > r*, completing the proof. □ 

We are now in position to prove our main theorem about the existence of invariant probability measures 
for NDDEs. 

Theorem 4.3. Let t be a positive constant, and B a nx n matrix with real entries and let g be in C^(M" x 
]R";R"). Consider the neutral delay differential equation given by 

^ {x{t) - Bx{t - t)) = g{x{t), x{t - r)), for t > 0, (4.73) 

on the phase space X :— C([—t, 0],R"). Assume that the dissipation conditions (4.65) and (4.66) are satisfied. 
Then, the continuous semigroup {S{t)}t>o acting on X which is generated by (4.64) possesses a global 
attractor A and the results in Theorems 2.1, 2.2 hold for {S'(t)}t>o- 

Proof. Keeping the previous notations, we note that for all £ X, = 9^'(0) — B(j>{—T) in this case. By 
Lemma 4.5, we infer the existence of an absorbing set for {S(t)}t>o, and since (4.64) is a particular form of 
the type of NFDEs handled by Lemma 4.4, to apply Theorem 4.2 and hence to infer all of the desired results 
we need now only check that D\ao < 0, where DIaq reduces simply here to sup {Re(A) : det I—Be^^'^ = O}. 

It is well known that in such a case Dq is stable (independently of the delay r) if and only if p[B) < 1 where 
p{B) is the spectral radius of B. This is the so-called Schur-Cohn condition; see [NicOl]. This last condition is 
guaranteed by the imposed dissipation condition (4.66) since ||i3|| < 1 necessarily, and p{B) < \\B\\ trivially. 
The proof is therefore complete. □ 

Remark 4.8. Note that this theorem may be extended to NDDEs of the form ^D^xt = G{xt) with Do as 
given in (4.59) and G being the functional representation associated to g in (4.64). The dissipation estimates 
and the stability criteria are however more involved. We leave such a generalization of Theorem 4.3 to the 
interested reader. 
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4.2.4 Application to a nonlinear Brayton-Miranker-like model 

We now return to a particular system of NDDEs, (4.76), arising from the transmission line problem discussed 
above in Section 4.2.1. Our goal is to provide conditions under which (4.76) exhibits a global attractor and 
hence to infer the existence of invariant measures for this example. In view of Theorem 4.3 we are left with 
the verification of the dissipation conditions (4.65) and (4.66) for this system which we establish below. 

Let 0<m, (j'<l;pSR, and let b and c be strictly positive real numbers. Let us introduce the notations 
$(t) := (01 (<), 02(0)^7 ^iid V = {vi,V2) and u = {ui,U2)- We use here bold typeface to distinguish vectors 
in M? from real numbers. We warn the reader that the symbol | • | will be used both for the absolute value 
of real numbers and for the Euclidean vector norm in R^. 

Let us consider Fi, ^2 : x — )• M two functions which satisfy respectively that there exist a' > 0, 
71 > 0, M[ > and Ml > such that for all (u, v) e ^ 



Mii^i(u, v) < —au\ + 7i, 
|?;2i^i(u,v)| < Miui+M[v2, 

and that there exist 72 > 0, M'2 > and AI2 > such that for all (u, v) G 



(4.74) 



M2-F2(u, v) < -a'ul + 72, 
|wiF2(u,v)| < M2W2 + M^wi. 

For instance Fi{u,v) = ~aiUi{l + |vp)^^ verifies the conditions (4.74) for i — 1 and (4.75) for i — 2. 
We consider now the following system of NDDEs: 

dt\ 



(4.75) 



'Mt) 




'0 q 




'Mt-r) 




P 




Miit) 




"Fi($(t),$(t- 


r)) 


02 (i). 




m 




>2(t-T)_ 











+ 


F2 ($(<), $(t- 





(4.76) 



We introduce furthermore k := max(m, q), which is the norm of the matrix arising in the LHS of (4.76) 
which has to be strictly less than 1 in view of Lemma 4.6, which explains the constraints imposed on m and 

q- 

The condition (4.65) to satisfy here can be written as: 



N{u,v) := {ui - qv2)(p - bui + Fi{u,v)j + (^2 - mwi)(^ - CU2 + F2(u, v)^ < 7 - a|up + ;3|vp, 

with a > 0, /? e M and 7 > to find independently of u, v e M^. 
An easy computation shows that, 

N{u, w) < 71 + 72 + pui + \pqv2 \ — hu\ — cu\ — a'|u|^ + A/iUi + M[vi + M2U2 + M2V2 + cmviU2 + bquiV2, 

now by applying the Young inequality on the two last terms and the e- Young inequality appropriately on 
the rest of the terms concerned, it is easy to show that for any e > there exists 7c > such that, 



A(m, v) < 7e 



- min(6, c) + a' — e) \u\ 



— max(&, c) + e) 



which shows that (4.65) is satisfied with B 



q 
m 



and 5(u, v) = (p - bui + i^i(u, v), -CU2 +i^2(u,v))^ 



Let us introduce i min(6, c) + a' — e and /?£ :— ^ max(6, c) + e. Assume now that a' is such that 

max(6,c) < (^ min(6, c) + a'). Let e > be fixed sufficiently small such that 2/3^ < a^. Now if we assume 

furthermore that k < —1 + ^ 2(1 — ^), then from Lemma 4.6 we deduce that there exists t* , such that 

(4.66) is satisfied for t > t*. 

By using now Theorem 4.3 we have thus proved the following proposition which constitutes, to the best of 
the authors' knowledge, the first result regarding existence of invariant measures for Brayton-Miranker-like 
models such as (4.76). 
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Proposition 4.3. Assume that a' arising in (4.74) and (4.75) is such that max(6,c) < (iniin(fe, c) + a'). 
For some e > such that 2/3^ < a^, assume furthermore that the coefficients q and m of the matrix arising 
in (4.76) satisfy: 

max(g, m) < - 1 + ^2 ^1 - 

where a ^ := i min(6, c) + Q!' — e and jS^ :~ i niax(&, c) + e. Let t* = — ^ log (P(max(q, m) + 2)/P(niax(g, to))) 
with P{x) :~ —[{x — 1)^ + 2(^ — 1)]. Then for all t > t* the continuous semigroup generated by (4.76) 
possesses a global attractor A living in C([— r, 0], M^). This semigroup possesses furthermore invariant prob- 
ability measures whose support is contained in A and which satisfy the "weak ergodic" property (4.63). 



Appendix: invariant measures and the global attractor 

For the sake of completeness we recall briefly a proof of the classical fact that invariant probability measures 
must always have their support contained in the global attractor; see e.g. [BCFM95, Lemma 4.2] or [FMRTOl] 
given in the concrete case of the 2D Navier-Stokes equations. 

Lemma 4.7. Let {S{t)}t>o he a continuous semigroup defined on a metric space {X,d). Suppose that 
{S{t)}t>o possesses a global attractor A. Then any invariant borel probability measure m (relative to 
{S{t)}t>o) has its support contained in A so that, in particular, m(^) = 1. 

Proof. For 5 > form the sets As — {y ■ mixeAd{x,y) < 5}. In view of the basic continuity properties of 
measures, it is sufficient to show that, m{As) = 1, for every 5 > 0, since, evidently, A = HsyoAs- 

To this end fix (5 > 0. Since A is attracting, for every R > and each x G X we may select tj^ > 
such that S(tfi)Bfi{x) C As, where Bfi{x) is the ball of radius R around the point x. This implies, 
Br{x) C S{tii)^^ S{t]i)BB,{x) C S{tfi)^^As- Thus, due to the invariance of m, we have, for every R > 
that m{Bji{x)) < m{S{t]i)~^As) = m{As)- Since for every x the collection {B]i{x)}b.>o is a nested collection 
of sets with Uij>o-B/?(a;) = X we now infer m{As) = 1 by again invoking basic continuity properties of 
measures. The proof is therefore complete. □ 
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